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Abstract 

The double-dipole strength distribution in °Ca is calculated microscopically 
within a model space of Iplh - and 2p2h excitations. Anharmonic effects in 
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OO , the centroid energies of the + - and 2 + components are found to be small, 
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in agreement with experimental observation. Firm conclusions about the 



t^j- | spreading width cannot be drawn, as yet, due to computational limitations in 

Q\ ■ 

^ I the number of 2p2h states. 



1. INTRODUCTION 



^ ■ In pion double-charge-exchange reactions as well as in peripheral heavy-ion collisions 



at relativistic bombarding energies > 1 GeV/nucleon the double-giant- dip ole resonance 
(DGDR) is strongly excited (see ]iy2[ for recent reviews). In the latter case there is a strong 
longitudinal focusing of the electromagnetic field in the target rest frame. This greatly en- 
hances the field intensity in the vicinity of the target nucleus thus increasing the probability 
for two-photon absorption from the ground state. The situation is somewhat similar to 
multi-photon excitations of atoms with intense laser pulses. 

Most simply, the DGDR is understood as a two-phonon excitation (with J n = + ,2 + 
in a spherical nucleus) which, in the harmonic limit, is located at twice the giant dipole 
resonance (GDR) energy. This expectation is largely confirmed by experiment [[J. On the 
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other hand, the single-dipole state 
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£>> = £mf|o>> 
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as well as the double-dipole state 



z z 



m> = (E^)(E^)I°> 



(2) 



i=i j=i 



are not eigenstates of the nuclear Hamiltonian and hence acquire a width. Here the exper- 
imental situation is less clear and values for Tdgdr are observed which are bracketed by 
\^2Tgdr an d 2Tgdr @- The lower value is expected from a folding of two Gaussians while 
the larger width results from the folding of two Lorentzians (recall that the Lorentzian gives 
a good fit to the total photoabsorption cross section). 

To reach a quantitative understanding of the DGDR characteristics a microscopic de- 
scription is needed. A first calculation within the second RPA formalism |||4| was presented 
in ref. [||] using a separable residual interaction and various degrees of approximation. More 
recently ||, particle- vibration coupling calculations of the semiclassical Coulomb excitation 
cross section have been reported which, by construction, select a restricted set of 2p2h dia- 
grams and have difficulties with properly imposing the Pauli principle. In the present paper 
we wish to present a calculation of the single- and double-dipole response function which 
diagonalizes a given model Hamiltonian, H, in the space of Iplh - and 2p2h excitations 
avoiding diagrammatic selections and fully respecting the Pauli principle. The theory has 
been used previously to provide a realistic description of linear response functions over a 
wide range of excitation energies ||. 

2. THEORETICAL DEVELOPMENT 

We start by expanding the wave packets \D) and \DD) (eqs. (1) and (2)) in exact 
eigenstates \N) of the nuclear Hamiltonian H as 



D)= 52(N\£>\0)\N) 
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\DD) = 52(N\DD\0)\N) (3) 
where D = Ylf=i ^% T i ■ When writing H as a mean field part and a residual interaction 

H = Y1 e i a \ a i + \j2 v ijM a\a)aia k . (4) 

i ij,kl 

( v ij,H denote antisymmetrized two-body matrix elements) the exact eigenstates \N) are 
obtained by diagonalization in the space of npnh excitations. A minimal truncation of this, 
prohibitively large, space which can accommodate both single- and double dipole excitations 
must clearly encompass the lplh and 2p2h sector, so that 

|ao = E^|i> + E4I2>- (5) 

1 2 

Due to the one-body nature of the dipole operator, the overlap matrix elements (iV|D|0) 
then project onto the lplh subspace, |1), with 

(N\D\0)=J2c»*{l\D\0) (6) 
i 

while (N\DD\0) project onto the 2p2h space, |2), with 



(N\DD\0) = J2<$* (21^10)- (7) 

2 

Knowing the expansion coefficients, cf 2 an d the eigenvalues of the Hamiltonian matrix 
in these spaces one then readily computes the 'strength functions' 

S D (E) = J2\(N\D\0)\ 2 5(E-E N ) 

N 

S DD {E) = Y,\{N\DD\Q)\ 2 5{E - E N ). (8) 

N 

The angular- momentum-coupled expression for Sod is somewhat involved and will be given 
elsewhere 0. For comparison with experiment it is useful to define the centroids and 
variances of the strength distributions (8) in terms of their energy moments 

m D(DD) = J dEE k S D{DD) (E) (9) 

as 

E G dr = m x D jm Q D \ a GDR = y/m 2 D - (m D ) 2 /m° D (10) 
and similarly for the double-dipole excitation. 



3. RESULTS 



Although not measured, so far, we have chosen the 40 Ca nucleus since the dimensions 
of the model space can be kept under reasonable numerical control. For the mean field we 
have chosen the same Woods-Saxon potential as in ref. |fl. Continuum excitations have been 
treated by discretization through an expansion of the positive-energy single-particle states 
in a harmonic oscillator basis. The density- dependent zero-range interaction of ref. M has 
also been adopted. Its form ensures the Pauli principle. To get a quantitative description 
of the GDR we have increased the isovector part, to 300 MeV fm 3 for, both, the interior 
and exterior part of Vqi. When truncating to Iplh and 2p2h configurations with excitation 
energies < 30 MeV (which yields 375 states) one obtains the strength distribution displayed 
in the upper part of Fig. 1. The centroid energy, E GDR is at 19.7 MeV as compared to 20.3 
MeV, deduced the global A-dependence E GDR = 31. 2A^ 3 + 20.6A- 1 / 6 |0|. Upon further 
increase of the model space, both the total transition strength and the centroid energy 
remain unchanged. For further comparison, the lower part of Fig. 1 shows the normalized 
photo absorption cross section || 

a,(E) = J2\(N\D\0)\ 2 E N 5(E - E N )/m l D (11) 

TV 

together with data from Compton scattering |TT| (for the theoretical distribution a 
Lorentzian smoothing of 1 MeV has been used). While not reproduced in all detail the 
overall agreement with experiment is quite satisfactory. 

In the harmonic limit, the double-dipole state is expected at 2 x Eqdr and a truncation 
at 30 MeV excitation energy is clearly insufficient. On the other hand the density of 2p2h 
states increases rapidly with excitation energy. A reasonable compromise is to truncate at 
45 MeV (where the number of states for the + -component of the DGDR is 1067 while 
for the 2 + -component one obtains 4144 states). At 45 MeV the total transition strength 
as well as the centroid energy are saturated. Any further increase in truncation energy 
will thus only influence the variance, (Jdgdr, i- e - the width. To test harmonicity one can 
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compare the DGDR centroid energies E^ GDR = 38.67 MeV and Ej^ GDR = 38.38 MeV to 
2 x Eqdr — 39.34 MeV, indicating that the anharmonicity is quite small. This is corrobo- 
rated by the fact that the + - and 2 + components are nearly degenerate. The normalized 
DGDR strength distributions, (Fig. 2), indicate some differences between the + - and 2 + 
components, however, in fine structure as well as in the overall width (the + -component is 
not strongly excited in experiment [|T2]j). 

To decide whether T DGDR /T GDR is closer to v2 or 2 we have determined the ratio of 
variances ugdr/^dgdr yielding 1.61 for the + -component and 1.24 for 2 + component (ctgdr 
was evaluated in a window from 10-30 MeV while (J^gdr was obtained in a window from 
30-50 MeV). It would be premature to conclude that a ratio of \[2 is favored, however, 
since the model space might still be too small to get a reliable answer of the width. We 
estimate that this limitation can be overcome. Before embarking on large-scale calculations, 
however, one should evaluate the Coulomb excitation cross section for direct comparisons 



with experiment [O] 



4. SUMMARY 

To get a quantitative assessment of the DGDR , which has been observed in several nuclei 
0, we have presented a microscopic calculation of the transition strength distribution in 
40 Ca within the space of lplh - and 2p2h excitations. By comparing the mean energies 
of the single- and double-dipole response it is found that anharmonicity effects are quite 
small, in good agreement with experiment. For the respective widths the conclusions are 
still hampered by computational restrictions of model space. Before enlarging the space it 
would be desirable to have a proper description of the Coulomb excitation cross section || . 
As another improvement, RPA ground-state correlations should be included. While this is 
trivial in the lplh sector, the 2p2h sector is more challenging numerically fi~3 |. 
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Figure Captions 

Fig. 1: upper part: the calculated dipole transition strength distribution in 40 Ca (eq. (8)). The 
high-energy tail has been multiplied by a factor of 10; lower part: the photoabsorption 
cross section, normalized to the TRK sum rule. The data, indicated by the filled 
squares, are inferred from the Compton-scattering analysis of ref . [9] . 

Fig. 2: The transition-strength distribution for the + -component (dashed line) and 2 + - 
component (full line) the DGDR in 40 Ca. For comparison the single-dipole strength 
function is also given short-dashed line). In all cases the discrete spectra have been 
folded with a 1 MeV-width Lorentzian (the experimental resolution is typically 1-2 
MeV [1]. 
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